######### EXERCICIOS 9 - ANA CAROLINA PAVAN ########

## QUE FRIO! ###
Celsius = function(x)
	{
	C = 5/9 * (x- 32)
	return(C)
	}
Celsius(74)

### Anlise exploratria simultnea de duas variveis ###
require(plotrix)
# ??multihist no foi muito elucidativo...
exploratoria= function (x,y)
	{
  	par(mfrow=c(2,2))
	multhist(list(x,y),main="Histograma de X e Y")
  	boxplot(x,y, main="Boxplot X e Y")
  	qqnorm(x,pch=16)
	qqline(x)	
	par(new=TRUE)  	
	qqnorm(y, ,pch=16, col=2, ann=FALSE, axes=FALSE)
	qqline(y,col=2)  	
	axis(side=4, at=y)
 	par(new=FALSE)
  	plot(y~x, main=paste("y em funo de x"))
  	par(mfrow=c(1,1))
  	return(list("Sumrio x"=summary(x),"Sumrio y"=summary(y), "Correlao x/y"=cor(x,y)))
	}
#Testando a funo:
x=rnorm(10,10,1)
y=rt(10,10,1.5)
exploratoria(x,y)

### ndices de Diversidade de Espcies ###

#ndice de Shannon: H=-?pilnpi
#ndice de Simpson: D=?pi2

ind.divers = function(x,indice="H")
	{
	abund.sp=apply(x,2,FUN=sum,na.rm=TRUE)
	pi=t(x)/abund.sp
	pi2=pi^2
	log=pi*log(pi)
	if(indice=="D")
		{
		#Diversidade de Simpson		
		valor=apply(pi2,1,FUN=sum,na.rm=TRUE)
		}
	if(indice=="H")
		{
		#Diversidade de Shanon
		valor=(-1)*apply(log,1,FUN=sum,na.rm=TRUE)
		}
	return(valor)
	}
# Validando a funo com dados meus:
dados<-read.table("dados Tocantinzinho.csv", header=TRUE, sep=";")
x<-as.matrix(dados)
ind.divers(x)
ind.divers(x,indice="D")

### FIM! ###






  